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Abstract. - We show that an electronic phase separation (EPS) transition described by the 
Cahn-Hilliard theory yields regions of low free energy forming grains of low and high charge 
densities. These local differences in the potential energy are studied numerically and used here for 
the first time as the origin of an attractive interaction that gives rise to local pair formation. The 
resistivity transition Tc occurs due to Josephson coupling among these superconducting regions or 
grains. Using this approach within the Bogoliubov-deGennes calculations we derive the Gaussian 
shape of Tc against the entire doping level p. We find that every pairing gap develops locally at a 
temperature Tp, following the relation 2A/kBTp ~ 8.0 in close agreement recent measurements. 



Introduction. — There are increasing evidences that 
the charge distribution in the Cu02 planes of the high 
temperature superconductors (HTSC) is microscopically 
inhomogeneous. Several different experiments like neu- 
tron diffraction [1-3], muon spin relaxation (fiSR) [4,5], 
NQR and NMR [6,7] have detected a position dependent 
electronic density. These experiments indicated that such 
disorder is stronger on the underdoped side of the phase 
diagram and it is possibly related with the non Fermi liq- 
uid behavior of the normal phase. However, recent STM 
studies on Bi2212 reveal spatial variations of the electronic 
gap amplitude on a nanometer length scale even on over- 
doped compounds [8-10]. These data on overdoped sam- 
ples, which behave as normal Fermi liquids at high temper- 
atures, indicates that a phase separation transition may 
occur at temperatures not much larger than Tc- A tran- 
sition at such low temperatures (of the order of lOOK) is 
likely to be an intrinsic electronic phase separation and 
not due to ionic mobility. 

On the theoretical side there are many predictions that 
hole carriers may segregate into hole-rich and hole-poor 
regions at low doping due to strong carrier-carrier corre- 
lation effect [11-13]. 

Here we deal with a completely different approach, 
namely, a phase separation transition driven by the mini- 
mization of the free energy [14]. As the temperature goes 
below the pseudogap temperature [15, 16] T*, the phase 



separation (PS) process starts. Such second order phase 
transition is a direct explanation to the NQR measure- 
ments [6] of two different signals coming from two differ- 
ent types of local doping domains as the temperature goes 
down. The origin of this EPS transition is the proxim- 
ity to the insulator AF phase, common to all cuprates, as 
we derived from the principle of the competing minimum 
free energy [14]. When the temperature decreases below 
Tps{p), the free energy of the homogeneous system with 
an average doping level or charge density p becomes higher 
than the disordered one, made mainly of two values of the 
local charge concentration at each point . According the 
stripe phase measurements [1,2] and the NQR data [6], the 
local values of the density p(ri) = p{i) follows a bimodal 
distribution [17] formed of AF domains with p{i) « and 
high hole density domains with p{i) « 2p. These domains, 
clusters or grains in the Cu-0 planes are of nanometer size 
containing 10-100 sites. 

The Electronic Phase Separation. To trace the 
EPS and the cluster formation we use the general theory of 
Cahn-Hilliard (CH) [18]. It describes how a system evolves 
from small fluctuations around the average charge concen- 
tration p near the phase separation temperature Tps{p) 
to a complete separation into low and high density grains, 
passing by intermediates configurations as the tempera- 
ture decreases. The order parameter of such transition 
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is the difference between tfie temperature dependent lo- 
cal charge or doping concentration p(i,T) and the aver- 
age doping level p, i.e., u{i,T) = (jpli,T) — p)/p. The 
Ginzburg-Landau (GL) free energy functional in terms of 
u{i,T) for a given compound p near the transition tem- 
perature is given by 



fii,T) = ^e^\Vuit,T)\^ + V{u{i,T)) = K + Vgl. (1) 

Where the potential VGL{p,i,T) = A^{T)u'^ /4,+ 
A^{T) = a{Tps{p) — T), a and B are temperature 
independent parameters, e gives the size of the grain 
boundaries among two distinct phases [19,20] . Vgl (p, «, T) 
changes with the site position "i"; V{u = 0) = at 
the grain boundaries where u — 0, and has two minima 
VGL{u,t) = —A'^{T)/B in the lowest and highest den- 
sity sites where the order parameter assumes the value 
Umin{T) — A{T)/ B. Thcsc two Vciminima regions 
work as attractors for the holes, giving rise to the clus- 
ters and to the main point here; it creates an effective 
two-body attractive potential. 

The CH equation can be written [21] in the form of a 
continuity equation of the local free energy /(i,T), dtu — 
—V.J, with the current J = MV{6f/6u), where M is the 
mobility or the charge transport coefficient. Therefore, 



dt ^ 



A^{T)u-B^u^ 



(2) 



We have already made a detailed study of the CH dif- 
ferential equation by finite difference methods [19] which 
yields the density profile u{p, i, T) in a 105 x 105 array as 
function of the time steps, up to the stabilization of the lo- 
cal densities, using parameters in the CH simulation that 
yield stripe [22,23] and patchwork [20,24,25] patterns at 
intermediate time regimes. 

The Fundamental Interaction. Here we intro- 
duce a new approach in order to derive the attractive po- 
tential that segregates the holes into grains of low and high 
density and which will be used as an effective two-body at- 
traction to calculate the superconducting properties. The 
justification to this procedure is because, given the meso- 
scopic size of a grain and the low doping values of the 
HTSC, there are very few holes in each confined region. 

Consequently we follow numerically the "kinetic" and 
potential energy (Vgl) map, that are the first and second 
term in Eq.(l). These energy maps show that the grains, 
either the low and high density ones, are regions of free en- 
ergy minimum, as shown in Fig.(l). This figure describes a 
common intermediate regime of disorder, between a homo- 
geneous system with a Gaussian distribution of densities 
and a bimodal distribution for a complete phase separa- 
tion, as it is shown in the inset by the histogram of local 
densities p{i). 

Following the Vgl values, from low temperatures up 
to Tps{p) when the grains disappear completely, we can 




Fig. 1: (color online) The Vgl potential energy map showing 
the regions of minimum potential separated by the lines of grain 
boundaries. The inset is the local charge density histogram 
for p = 0.16 showing the tendency toward a bimodal charge 
distribution. The points marked below the inset are where we 
made the local density of states calculations presented in the 
end. 



obtain the qualitative behavior of this attractive poten- 
tial. Here we use parameters to reproduced the Bi2212 
low temperature data of McElroy et al [8] 

VGL{p,hT) = V{p)V{T) = (-0.9 + 2.8p)|«(z,r)p, (3) 

where the values are in eV , V{p) is linear and vanishes at 
p « 0.32 following the behavior T*(p) or Tpsip). V{T) 
vanishes at Tps{p) and increases as T decreases. The 
u{i,T) is also parameterized in a similar way, |m(i,T)|^ = 
u{i){\ - r/Tps)(3-i'/i'ps) showing how the system be- 
comes homogeneous at Tps, and the other limit at T = 
of largest phase separation with A = B = 1. 

The energy minimum at the grains and the grain bound- 
ary energy barrier keep the holes confined forming metallic 
clusters with typical Fermi energies much smaller than the 
bulk Debye energy, the so called anti-adiabatic regime [26] . 
This condition with the effective hole attraction toward 
the cluster center is highly favorable to bipolaron forma- 
tion [26]. Then the effective hole attraction, modelled by 
Vgl, is taken as the origin of the superconducting inter- 
action in the form of a local two-body attraction in the 
Bogoliubov-deGennes (BdG) theory. This assumption is 
in agreement with the observation that a large electron 
mobility and a small Fermi energy produces a large Nernst 
Effect [27] as observed in many cuprates [28]. 

BdG-CH Combined Calculations. To calcu- 
late the local superconductivity amplitude or gap func- 
tion Ad{i,T) in a plane, we use the BdG theory with 
the extended Hubbard Hamiltonian in a similar fashion 
as we did before for a phenomcnological next neighbor 
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potential [20,22-25]. Except from the temperature de- 
pendent Vgl{p- 'i-,T), introduced here for the first time, 
and the granular charge density taken from the CH cal- 
culations and used as an input, all the others parame- 
ters are similar to values previously used in other calcula- 
tions [14,20,22,23,25]. 

The BdG-CH combined calculations yield large super- 
conducting amplitudes in high density grains and low or 
almost zero local amplitudes at regions with low densities. 
We define the local superconducting temperature Tc{i) as 
the onset temperature for Ad{i,T). The largest value of 
Tc{i) in a given compound determines the temperature 
Ton{p) which marks the onset of superconductivity. Since 
Ton{p) is correlated with the potential V{p,i, T) it is larger 
in the underdoped region and decreases almost linearly as 
the average doping increases, similar to the Nernst effect 
onset temperature [22,28] or a signal above Tc associated 
also with pair formation by tunneling measurements [29]. 

As the temperature decreases below Ton{p) some grains 
become superconductors, the zero resistivity transition 
takes place when the Josephson coupling Ej among these 
grains is sufficiently large to overcome thermal fluctua- 
tions [30]. Thus Ej{p,T = Tc) « ksTcip) what leads 
to phase locking and long range phase coherence. Conse- 
quently the superconducting transition in cuprates occurs 
in two steps, similar to a superconducting material em- 
bedded in a non superconducting matrix [30]: First by the 
appearing of intragrain superconductivity and by Joseph- 
son coupling with phase locking at a lower temperature, 
what provides a clear interpretation to the presence of the 
superconducting amplitude Ad{i,T) above Tc{p). 

The Cuprates Phase Diagram. — In this approach 
Tc{p) is not directly related with the intragrain super- 
conductivity, and the amplitudes Ad{i,T) do not change 
appreciably around Tc{p), specially for underdoped com- 
pounds, characterized by very large To„(p). This conse- 
quence of the present granular theory has been experimen- 
tally verified by temperature dependent tunneling [31] and 
angle resolved photon emission [32]. Using to the theory 
of granular superconductors [33] to the electronic grains. 



Josephson Coupling and kgT 



Ej{p,T) = 



nh 
Ae^Rn 



(4) 



Where A{T,p) is the average of the BdG superconducting 
gaps calculations Ad{i, T) on a 24 x 24 square taken from 
the 105 X 105 mesh after the CH simulations as that shown 
in Fig.(l). The Rn is the normal resistance of a given 
compound, which we take as proportional to the pab mea- 
surements [34] on the complete series of La2-pSrpCu02- 
These Rn values are given in the legend of the Fig. (2) top 
panel. The average A{T,p) as function of p are plotted in 
the low panel of Fig. (2) with the Tc{p) results from Eq.4) 
are in the inset and yields the well known dome shape with 
excellent agreement with the Bi2212 values. This is one of 
the most important result of our CH-BdG calculations. 
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Fig. 2: In the top panel we plot ksT and the Josephson cou- 
pling among superconducting grains E,j(p, T) for some selected 
doping values as function of T. The curves intersections give 
the dome shaped Tc{p), as plotted in the inset in the down 
panel where we show the average A(T,p) calculations. 



The STM Results and Interpretation. Accord- 
ing to the recent STM data [9, 10], the low values of 
T*(p) ~ Tpsip) for overdoped samples rules out any ionic 
mobility as the origin of the inhomogcneitics in overdoped 
samples and probably at all dopings. Thus we want to 
show that these STM results can be interpreted by the 
granular behavior resulting from the EPS, calculating the 
symmetric local density of states (LDOS) at difi^erent local 
doping, 

Ni{E)=Y,[ I «„(x,)|2 + |t;„(x,)n X 

n 

[fUE - En) + fi{E + En)], (5) 

where /„ is the Fermi function, the prime is the deriva- 
tive with respect to the argument, and Un,Vn and En 
are respectively the eigenvectors and positive eigenvalues 
(quasi-particles exciting energy) of the BdG matrix equa- 
tion [14,20,22,23,25]. 

Here we concentrate on a p = 0.23 compound that is 
close to the Bi2212 compounds used in recently STM ex- 
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Fig. 3: In the panel a) we plot the LDOS at a "metallic" 
site (i=175) inside a high density grain with A(r — 30K) — 
17meV. As the temperature increases the coherent gaps close 
(at T « 52K). In panel b) we show a LDOS with at an "insu- 
lator" grain (i=109) which has both p{i) and A(r = 0) almost 
zero. At T = 120K both LDOS converges to the same value 
all over the plane because the proximity to Tps ~ 140K. The 
inset shows the low temperature STM data of McElroy et al [8] 
displaying the same behavior of our calculations. 



periments with many novel results [9, 10]. The low tem- 
perature CH EPS calculations on this sample with the es- 
timated value of Tp5 = 140K yield that « 20% of the sites 
are in very low densities regions with p(i) < 0.05. These 
low density unitary cells have quite distinct LDOS behav- 
ior than the high density ones as it is shown in Fig. (3). 

The calculations at sites in the high density regions have 
larger superconducting amplitudes A{i,T) with well de- 
fined (coherent) and high LDOS peaks. A(i,T) is cal- 
culated directly from the BdG equations but can be also 
estimated by the position of the first LDOS peak. As the 
temperature increases, it shows the closing of these peaks 
and the building up of spectral weight at the Fermi level 
(E=0) (Fig. (3a)). In particular, for this case, we have a 
A(0, i — 175) « 17.7meV and the gap derived from the 
peaks decreases and closes near T = 52K. This behavior 
is quite common to the sites with p{i) > 0.05 and with 



Tc{i) around 45K and 60K. 

On the other hand, the calculations at points with 
p{i) < 0.05 display a rather different behavior as it is 
shown in panel b). As an example we show the LDOS at 
i = 109 with p{i) « and A(0, i = 109) « 0. It shows a 
much smaller LDOS (compare the vertical scale of panel 
a) and b)) without the low temperatures coherent peaks as 
in panel a) and with small oscillations far from the Fermi 
energy given an impression of a larger (incoherent) gap. 
These small oscillations followed by other peaks are quite 
similar to the firsts few peaks in the work of McElroy et 
al [8] showing in the inset of panel b) . As the temperature 
increases and reaches T — 5QK we see that the LDOS at 
panel a) is almost close yielding & A{i — 173,t — 5QK) ^ 
while the LDOS at panel b) has still a large dip that gives 
a A{i = 109, t = 50i^) ~ 20meV. This dip in the LDOS 
stays much above the resistivity transition Tc = 60K in 
close agreement with the STM maps of Gomes et al [9] . 

Since the Vol potential (Eq.(3)) vanishes with the or- 
der parameter, the LDOS for both cases converge to the 
same value as the temperature approaches the "melting" 
temperature Tps — 140K, as it is demonstrated in the 
Fig. (3) by the T = 120K curve. 

Another striking result of the CH-BdG calculations is 
that, despite the uncertainty on Tc{i) (or Tp in the nota- 
tion of Gomes et al [9]) for very small gaps, mostly of our 
results follow close the measured relation 2Ad/ KBTc{i) « 
7.9 [9]. For the case of i=175, Ta{i) = 52K = 4.5meV 
and 2Ad(i,0) = 34mev which gives a ratio of 7.6. 
For i=109, taking the first peak in Fig. (3b) at 26meV 
which closes completely at T w 74K, we obtain again 
2A/ KBTc{i) ~ 8.1 again in close agreement with the STM 
data [9] but larger than the value measured by tunneling 
of 2Ad/KBTc{i) « 6.0 [29]. Since the lower density (insu- 
lator) grains have much smaller LDOS around the Fermi 
level, they also have lower local conductivity, although the 
gaps measured by the position of the LDOS peaks, as dis- 
cussed above, are larger than high density sites with their 
coherent peaks. Consequently regions with larger gaps 
have lower conductivity as it was recently measured [10]. 

Conclusion. — In summary we have made a detailed 
study of the EPS in HTSC using the CH theory that 
allows us to followed the local free energy minima that 
generate the high and low density grains in the Cu-0 
planes. The differences in these free energy local min- 
ima were calculated numerically and used in the BdG 
approach as the superconducting interaction in the elec- 
tronic grains of cuprates in the non adiabatic limit. These 
calculations give rise to the intragrain superconductivity 
and they provide a scenario to the pseudogap phase, as 
composed of local regions with finite superconducting am- 
plitude A{i,p,T) without phase locking. The Josephson 
coupling calculations in connection with the measured val- 
ues of the resistivity yield an accurate Tc{p) curve for the 
whole doping values of the Bi2212 series. The LDOS cal- 
culated in the high density grains yield the coherent peaks. 
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and those in the low density regions give larger gaps and ill 
defined peaks in agreement with the STM data in Bi22f 2. 
As far as we know, this present work is the only one to 
give an interpretation to all the whole STM results. 

All of these calculations in close agreement with current 
data led us to conclude that the EPS is an important 
ingredient above the superconducting phase and generates 
mostly of all the intricate normal phase physics of the 
HTSC. 
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[31] Suzuki Minoru and Watanbe Takao , Phys. Rev. Lett, 

85 (2000) 4787. 
[32] Kanigel A. ET al, Phys. Rev. Lett, 101 (2008) 137002. 
[33] Ambeogakar V. and Baratoff A., Phys. Rev. Lett. , 

10 (1963) 486. 

[34] Takagi H. et al, Phys. Rev. Lett. , 69 (1992) 2975. 



REFERENCES 

[I] Tranquada J.M. et al, Nature (London), 375 (1995) 
561. 

[2] BiANCONi A. ET AL, Phys. Rev. Lett, 76 (1996) 3412. 
[3] BoziN E.S. ET AL, Phys. Rev. Lett, 84 (2000) 5856. 
[4] Uemura Y.J., Sol. Stat Comm. , 126 (2003) 23. 
[5] SONIER J.E. ET AL, Phys. Rev. Lett , 101 (2008) 117001. 
[6] Singer P.M., Hunt A.W., and Imai T, Phys. Rev. Lett 

, 88 (2002) 47602. 
[7] Grafe H.J. ET AL, Phys. Rev. Lett, 96 (2006) 017002. 
[8] McElroy K. , ET AL, Phys. Rev. Lett, 94 (2005) 197005. 
[9] Gomes Kenjiro K. et al. Nature , 447 (2007) 569. 
[10] Pasupathy Abhay N. et al. Science , 320 (196) 2008. 

[II] GORKOV L.P. and Sokol A.V., JETP Lett , 46 (1987) 
420. 

[12] YuKALOV V.I. and Yukalova E.P., Phys. Rev. B, 70 
(2004) 224516. 

[13] Emery V.J. and Kivelson S.A., Physica C, 209 (1993) 
597. 

[14] de Mello E.V.L, Passos C.A.C and Kasal R.B., J. 

Phys.: Condens. Matter, 21 (2009) 235701. 
[15] TiMUSK T. and Statt B., Rep. Prog. Phys. , 62 (1999) 

61. 

[16] Tallon J.L. and Loram J.W, Physica C, 349 (2001) 53. 
[17] DE Mello E.V.L. et al, Phys. Rev. B, 67 (2003) 024502. 
[18] Cahn J.W. AND Billiard J.E. , J. Chem. Phys, 28 
(1958) 258. 

[19] de Mello E.V.L. , and Silveira Filho Otton T. , 
Physica A, 347 (2005) 429. 

[20] DE Mello E.V.L. and Caixeiro E.S., Phys. Rev. B, 70 
(2004) 224517. 

[21] Bray A.J. , Adv. Phys., 43 (1994) 347. 

[22] de Mello E. V. L. and DiAS D. N., J. Phys. CM. , 19 
(2007) 086218. 

[23] DiAS D. N. et al, Physica C, 468 (2008) 480. 

[24] DE Mello E.V.L. et al. Proceedings of SCES08, Phys- 
ica B (2009) . 

[25] Caixeiro E.S, de Mello E.V.L., and Troper A,, 

Physica C, 459 (2007) 37. 
[26] de Mello E.V.L. and Ranninger J., Phys. Rev.B, 55 

(,) 1997 14872 

[27] Behnia Kamran, J. Phys. -.Condens. Matter, 21 (2009) 
113101. 

[28] Wang Yayu , Li Lu and Ong N. P., Phys. Rev. B , 73 
(2006) 024510. 

[29] Mourachkine a.. Mod. Phys. Lett. B, 19 (2009) 743. 
[30] Merchant L. et al, Phys. Rev.B, 63 (2001) 134508. 



p-5 



